Exact solitary wave solutions for a discrete X<j)^ field theory in 1 + 1 dimensions 
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We have found exact, periodic, time-dependent solitary wave solutions of a discrete (f)^ field theory 
model. For finite lattices, depending on whether one is considering a repulsive or attractive case, 
the solutions are either Jacobi elliptic functions sn(a:, m) (which reduce to the kink function tanh(a;) 
for m 1), or they are dn(a:,m) and cn(2;,m) (which reduce to the pulse function sech(a;) for 
m — > 1). We have studied the stability of these solutions numerically, and we find that our solutions 
are linearly stable in most cases. We show that this model is a Hamiltonian system, and that the 
effective Peierls-Nabarro barrier due to discreteness is zero not only for the two localized modes but 
even for all three periodic solutions. We also present results of numerical simulations of scattering 
of kink-anti-kink and pulse-anti-pulse solitary wave solutions. 



PACS numbers: 61.25.Hq, 64.60.Cn, 64.75.-|-g 



I. INTRODUCTION 



Discrete nonlinear equations arc ubiquitous and pj^ 
an important role in diverse physical contexts P, 0. 
Some examples of integrable discrete equations include 
Ablowitz-Ladik (AL) lattice |3| and Toda lattice 
Certain non-Abelian discrete integrable models are also 
known 0. There are many non-integrable discrete non- 
linear equations such as discrete sine-Gordon (DSG) la, 
0], discrete nonlinear Schrodinger (DNLS) equation Q 
and the Fermi-Pasta-Ulam (FPU) problem Q. DSG 
is a physical realization of the dynamics of dislocations 
in crystals where it is known as the Frenkel-Kontorova 
model ^3- It also arises in the context of ferromag- 
nets with planar anisotropy [Tll | , adsorption on a crystal 
lattice and pinned charge-density waves jl2l |. Similarly, 
DNLS plays a role in the propagation of electromag- 
netic waves in doped glass fibers [13 and other optical 
waveguides and it describes Bose-Einstein conden- 
sates in optical lattices FPU has served as a fertile 
paradigm for understanding solitons, discrete breathers, 
intrinsic localized modes, chaos, anomalous transport in 
low-dimensional systems and the fundamentals of statis- 
tical mechanics Q . A discrete double well or discrete (j)^ 
equation is a model for structural phase transitions [16|, 
and may be relevant for a better understanding of the 
collisions of relativistic kinks [T7| . 

Obtaining exact (soliton-like) solutions is always desir- 
able, particularly for discrete systems where the notion 
of a discreteness (or Peierls-Nabarro) barrier 0, is 
an important one in that its absence is a likely indica- 



tion of integrability, e.g. in the AL case In addition, 
exact solutions allow one to calculate certain important 
physical quantities analytically as well as they serve as 
diagnostics for simulations. Recently derived summation 
identities involving Jacobi elliptic functions [21I |23 | 
led to exact periodic solutions of a modified DNLS equa- 
tion |23l |. In this paper we exploit similar identities to 
obtain exact periodic solutions of the discrete cf>^ model. 

The paper is organized as follows. In the next section 
we summarize the exact solitary wave solutions of the 
continuum model. In Sec. IIIII by identifying the rel- 
evant elliptic function identities, we derive a discrete 
model which allows for exact solutions. We then obtain 
these solutions and study their stability. In Sec. IIVI we 
explicitly write down the corresponding Hamiltonian dy- 
namics using a modified Poisson bracket algebra. Sec. Ivl 
contains numerical results for the scattering of both kink- 
and pulse-like solitary waves. Kink-anti-kink collisions 
appear to create a breather with some radiation, and 
pulse-anti-pulse collisions lead to a flip and little radia- 
tion. In Sec. IVII we compute the energy of the solit ary 
waves, and show that the Peierls-Nabarro barrier [la, [13 
due to discreteness is zero not only for the two local- 
ized solutions but even for the three periodic solutions. 
Finally, we summarize our main findings in Sec. IVIII 



II. CONTINUUM SOLITARY WAVES 

The double- well potential with the coupling parameter 
A and the two minima at (h — ±a 



y = ^ (02 - ar 



(2.1) 
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leads to the following relativistic field equation 



A0(r - a^) = , 



(2.2) 



2 



For A > 0, if one assumes a moving periodic solution with 
velocity v in terms of the Jacobi elliptic function sn(x, m) 
with modulus m 



4>{x, t) = Aa sn [/?(a; + xq — vi), m] 
one finds that 



(2.3) 



(2.4) 



dt"^ dx^ 

= -f3^a (1 - w^) [(1 + m) sn - 2m sn^] , 

where we have suppressed the argument of sn. Match- 
ing terms proportional to sn and sn'^ we find we have a 
moving periodic solution with 



(1 - + m) 

2m 



1/2 



1 + m 



(2.5) 



The usual kink solitary wave is obtained in the limit 
TO ^ 1. Since 



sn(a;, 1) — tanh(a:;) 



(2.6) 



we obtain 
(f>{x, t) = atanh-^ 



Aa2 



ll/2 



2(1-1-2) 



{x + xo~ vt) ) . (2.7) 



If instead A < 0, then there are dn Jacobi elliptic function 
solutions to the field equations. Assuming 

(j){x, t) = Aa dn [[5{x + xq ~ ut), m] , (2.8) 

and matching terms proportional to dn and dn^ . we find 
a moving periodic pulse solution with 



A = 



-Aa2 



(1-w2)(2_to) 



1/2 



2 - m 



(2.9) 



In fact, there is another pulse solution in terms of the 
Jacobi elliptic function of cn type. Assuming 



{x, t) = Aa cn [/3(.t + — wt), to] 



(2.10) 



Matching terms proportional to cn and cn^ we find we 
have a moving periodic solution with 



A = 



-Aa2 



1/2 



(1- w2)(2to- 1)_ 
2to 



2to- 1 

Note that this solution is valid only if m > 1/2. 



(2.11) 



The usual pulse solitary wave is obtained in the limit 
m —^ 1. Since 



dn(a::, 1) = cn(a;, 1) = sech(a::) , 



(2.12) 



we obtain 

4>{x, t) — a-\/2sech 



-Aa2 



2 \ 1/2 



(x + xo — vt) 



(2.13) 

Note that similar solutions and their stability were con- 
sidered by Aubry in the context of structural phase tran- 
sitions [2J|. In the rest of this paper we will refer to the 
sn solutions as kink-like, and to the dn and cn solutions 
as pulse-like. 



III. DISCRETIZATION OF FIELD THEORY 



A naive discretization of the field equation above is 



1 



-0„-h^(<?!)„+i+0„_i-20„) + A^„(a^-(/),l) = 0, (3.1) 

where e is the lattice parameter and the overdots repre- 
sent time derivatives. However, this equation does not 
admit solutions of the form 

0„(<) = Aa sn(/3[(n + c)e - wt], to), (3.2) 

where c is an arbitrary constant. To find solutions, one 
has to modify the naive discretization. The key for un- 
derstanding how to modify the equation comes from the 
following identities |23| of the Jacobi elliptic functions 
sn: 

TO sn(a;, m)^ [sn{x + I3,m) + sn{x — /3, to)] (3-3) 
= ns^ (/3, m) [sn(x + I3,m) + sn{x — (3, to)] 
— 2cs(/3, TO)ds(/3, TO)sn(a;, to) , 

cn: 

TO cn(a;, m)^ [cn(x -|- /3, m) + C'D-{x — /?, to)] (3-4) 
= — ds^(/3, m) [cn{x + /?, to) 4- cn(a; — /?, m)] 
+ 2cs(/3, TO)ns(/3, m)cn(x, to) , 

and dn: 

dn(2;, to)^ [dn(a; -|- /3, to) + dn(a; — /3, to)] (3-5) 
= — cs^ (/3, to) [dn(x -f- 13, m) + dn(.T — /?, to)] 
-|- 2ds(/3, TO)ns(/3, TO)dn(a;, to) , 



where 



1 



sn(a;, to) 



cn(x,TO) dn(a;,m) 
cs = — -, r , ds = 



sn(x,TO) ' 



sn(a;, m) 



(3.6) 



3 



It will be useful later to rewrite these identities in the 
form: 

sn{x + P , m) + sn{x — P , m) (3-7) 
= 2 cs(/3,m) ds(/3,77i) sn(x,TO) 
X [ns^(/3,m) — m sn(a;, m)^]""'^, 



cn(a; + /3, to) + cn(.T — /3, m) 

= 2 cs(/3,to) ns(/?,77i) cn(x,TO) 
X [ds^(/3, to) + TO cn(a;, m)^]^"'^, 



(3.8) 



and 



dn(x + /3, m) + dn(a; — /3, to) (3-9) 
= 2 ds(/3,TO) ns(/3, to) dn(a;,TO) 
X [dn(a;,TO)^ +cs^(/3, to)]"\ 

For the sake of brevity, in what follows we will suppress 
the modulus to in the argument of the Jacobi elliptic 
functions, except when needed for added clarity. 

A. Static lattice solutions 

Consider the naive lattice equation: 

1 {(f>n+i + (j)n-i - 20„) + XMa^ - 0n) = . (3.10) 

Here = 0[/3e(n + c),to]. A general ansatz to modify 
the solution is to multiply the second difference operator 
by the factor (1 — ct<Pn): with a chosen so that we get a 
consistent set of equations. That is, we consider instead 
the equation: 



1 — ao 



fn+l 



+ (f>n-l-2cbn)+XMa^-^l) =0- (3.11) 



Using a to eliminate the (p^ term leads to the result 

(3.12) 

This implies, in the static case, that the lattice equation 
is just a smeared discretization of the 0^ term. Namely, 
one merely needs to study the lattice equation: 



— {4>n+l + 4>n-l - 20„) = ^^((/)„+l + (t)n-l) ~ A0„ 



(3.13) 



. . , 2 

Assume a solution of the form (for A > 0) 

0„ = Aa sn„ , (3.14) 



where sn„ denotes sn[/3e(n + c),to] with c being an ar- 
bitrary constant. Note that we only need to consider c 
between and 1/2 (half the lattice spacing). On match- 
ing the coefficients of sn„ and sn„+i -I- sn„_i we obtain 



2 2 _ (1 +™) sn^iPe) 



Ae2 



(3.15) 



and 

Xa^ = § [1 - dn(/3e)cn(^e)] 
In the continuum limit (e ^ 0), we get 
Xa^ ,9 2to 



f3' 



1 -|- TO 



1 + TO 



(3.16) 



(3.17) 



which agrees with our continuum result (2.5) with w = 0. 
Also as TO — > 1 we recover the usual kink solution 



(/)„ = a tanh[/3e(7i + c)] . 
If instead, we have A < and we assume 
(hn = Aa dn„, 



we obtain: 



and 



2„2 



A'a 



Xa' 



2m sn^(/3e) 
A?cn2(/3e) 

_ dn{(3e) 
~ cn2(/3e) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



In the limit when the lattice spacing e goes to zero, we 
obtain 



Xa^ 
' 2-m 



A' = 



2-to' 



(3.22) 



same as the continuum case (2.9) with v ~ 0. 

For A < we also have a possible solution of the form 



we obtain: 



and 



2„2 



A' a 



Xa^ 



Aa cn„, 

2to sn2(/3e) 
X^dii^Pe) 

, cn(/3e) ■ 



dn^(/3e) 



(3.23) 



(3.24) 



(3.25) 



Again taking the lattice spacing to zero, we obtain 
Xa^ . , 2to 



2to- 1 



A' = 



2m - 1 ■ 



(3.26) 



same as the continuum case (2.11) with v ~ 0. Note that 
as before this solution is only valid if to > 1/2. For both 
these solutions, as to — > 1, we get our previous continuum 
solution: 



V2a sech[pe{n + c)] . 



(3.27) 
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FIG. 1: (Color online) Kink-like (A > 0) case: N-dependence 



of the lowest eigenvalue, uji, for A 
values of the elliptic modulus m. 
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FIG. 2: (Color online) Kink-like (A > 0) case: N-dependence 
of the lattice spacing, e, for A = 1, a = 1, and various values 
of m. 



FIG. 3: (Color online) Pulse-like (A < 0) dn case: N- 
dependence of the lowest eigenvalue, ujf, for A = —1, a = 1.5, 
and various values of m. 
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FIG. 4: (Color online) Pulse-like (A < 0) dn case: N- 
dependence of the lattice spacing, e, for A = —1, a = 1.5, 
and various values of m. 



B. Stability of solutions 

Let us now discuss the stability of these static solu- 
tions. To that purpose, let us expand 

0n = (f>i^^ + VAi exp(-icji), 

(s) . 

where is the known solitary wave exact solution and 
ipn is a small perturbation. Then we find that to the 
lowest order in ipn the stability equation is 

-AV„0i^)[0iti + '^i'^-li] 

+AV'„a2- ^(V-n+i+Vn-i) = 0. (3.28) 



Using the identities (|3.7(l . (3.8) and H3.9() . the combina- 

(s) (s) 

tion -I- 4>n-i] can be written as a function of 0„ 

alone. Thus schematically we have 

{uj^-Q^„ + (7„(V'„+i+^„-i) = 0, (3.29) 

or 

^pn+l + Ipn-l + hni>n = , (3.30) 

with 



gn 

where /„ and gn are well defined functions. We also have 
that the solutions are periodic on the lattice with N sites 
(jpn+N = 'fpn)- If we write the system of N equations in 
matrix form, as 

A [^] = [^], (3.32) 



5 



then the condition for nontrivial solutions is that 



One easily verifies that 



Det |A 



= 0. 



(3.33) 



In our simulations we require that the solution has ex- 
actly one period in the model space. Therefore, the equa- 
tions for f3 and A are supplemented by the equations 

NPe = 4 K{m) , (3.34) 

for solutions, 4>n ~ sn(/3ne) or 0„ = cn(/3ne), and 



NPe = 2 K{m) , 



(3.35) 



when (pn = dn(/?7ie). Here K(m) is the complete ellip- 
tic integral of the first kind 0, Hence, once the 
parameters m, A and a are specified, we need to solve 
a system of equations for /3, e, and A, i.e. we have to 
solve Eqs. l|TT31IXT51l?3^ or Eqs. 12201 Em or 
Eqs. (|3.24l 13.251 13.34|l . for the solutions sn, dn, and cn, 
respectively, for f3, e, and A as a function of N. 

Typical results are depicted in Figs. ^ and |21 for kink- 
like sn solutions, and in Figs.|21and21 for pulse- like dn so- 
lutions, respectively. Similar results can also be obtained 
in the case of the cn-type pulse solution. For the case 
of the kink-likc solutions and our choice of parameters, 
we find from Fig. ^ that stability requires N > 7 lattice 
sites. In the pulse-like case, we find stability for arbitrary 
values of TV. The magnitude of the lowest eigenvalues in- 
creases with N in the kink-like case, while it decreases 
with N for pulse-like solutions. For fixed m. A, and A, 
the lattice spacing, e, is not an independent quantity, but 
it is a well-defined function of N. As seen from Figs. |21 
and 21 the lattice spacing, e, is always a decreasing func- 
tion of N. 



IV. HAMILTONIAN DYNAMICS 

In this section we demonstrate that our discrete model 
is a Hamiltonian system. The equation which the static 
solutions obey, Eq. (|3.11|l . can be written as: 

1 (0„+i + - 20„) + X^^-^ = 0. (4.1) 

We recognize that this equation can be derived from the 
potential function 



(0n+l - (j^ny 

2e2 



(/-„(a2 -4,1) 



1 — ao 



Performing the integral we obtain explicitly: 

i>n+l - (t)nY 



Vr 



2e2 



(4.2) 
(4.3) 



2a 



2a2 



■n+l + 4>n-l - 24>., 
2 j,2 



(4.4) 



The continuum limit a — > (or e — > 0) of the last two 
terms in Eq. H4.3|l is given by 



j{9n -a ) - -a 



(4.5) 



We will now show that if we want our class of solutions 
to be the static limit of a Hamiltonian dynamical system, 
that we will be able to have single solitary waves obey a 
simple equation, but general solutions will obey a more 
complicated dynamics with terms proportional to (j)'^. For 
simplicity we will assume that the Hamiltonian takes the 
form: 



H 



E 



(4.6) 



with V given by Vq in Eq. H4.3|l plus possibly some 
additional terms that vanish in the continuum limit. 
Here 7r„ is the conjugate momentum and g[(j)„] a weight 
function. For generality we will assume, as in the 
case of the discrete nonlinear Schrodinger (DNLS) equa- 
tion I2J [2^ , that an extended Poisson bracket structure 
exists, namely: 



{0m, 7r„} = Sn,nf[4'n 



(4.7) 



and 



0n = {4n,H} = {(/)„, 7r„J = /[0„]-^, 

T^n = {7r„,i?} = 77— {7r„,0m} = -f[4>n]-^- 



d(t)n 



(4.8) 



From our ansatz, Eq. 1)4. 6|l . we obtain the first order equa- 
tions: 



and 



4>n = 7r„/[(/)„]g[0„] = 7r„/l[0„ 
.r , J-^ldg , dV 



(4.9) 



■ (4-10) 



This leads to the following second order differential equa- 
tion for 4>n 



1 fU IhU l^^^i2^1 dh f dg 

4>n = -j[(Pn\h[(l)n\Tr2 



h d(j)n 2h d(j)r, 



(4.11) 



For this equation to have the previously found static soli- 
tary waves, as well as the correct continuum limit, we 
need only that 



/[0n]/i[0n] = f[(f>n]g[<t)n] = 1 - a^l 



(4.12) 
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The two simplest cases are: (a) / = 1 (ordinary Poisson 
brackets) and (b) h = 1 (extended Poisson brackets). 
First consider the case / = 1. This requires 



h[(f>n] = g[(t)n] = 1 - a(t)l 



(4.13) 



From this we get, if we choose V — Vq with Vq given by 
Eq. H4.3|l . the equation: 



1 



20n) 



(4.14) 



1 



For a = Ae^/2 this equation has the previously found 
static lattice solitary wave solutions as well as the correct 
continuum limit. For the case h ~ 1, this leads to 



(4.15) 



We find that for V = Vo the equation of motion again 
becomes the same as Eq. (4.14). 

Unfortunately the presence of the (jP term does not 
allow for single elliptic time dependent solitary wave so- 
lutions that are of the form sn[/3{x — ct)], cn[/3(x — ct)] 
or dn[/?(a; — ct)] where x = ne. This is because the sec- 
ond derivative of the sn and dn (or cn) functions contains 
both linear as well as cubic terms, and the quantity 0^ is a 
quartic polynomial in sn or dn (or cn) . Thus the last term 
is equivalent to a non-polynomial potential when applied 
to a single elliptic function solution. Therefore in order 
to obtain a simple elliptic function solution in the time 
dependent case, one must add a non-polynomial poten- 
tial to Vq which is chosen to exactly cancel the last term 
when evaluated for a single elliptic solitary wave solution 
of the form sn[P{x — ct)], cn[/3(a- — ct)] or dn[f3{x — ct)]. 
This means that to obtain an equation that is derivable 
from a Lagrangian or a Hamiltonian we should add to the 
static potential Vq an additional contribution AV such 
that 



(1 - 



dAV 

d(j)n 



1 - a02 



(4.16) 



where is a single solitary wave described by a time 
translated elliptic function. It will turn out that AV 
needed to obtain a simple solution of the elliptic kind 
will depend on the velocity squared of the solitary wave 
as well as the type of solution (pulse- or kink-like). 

In general we have that the Hamiltonian dynamics 
leads to 

Since fh = 1 ~ a(j)f^ and h = /g, it follows that 
1 dh f dg a(j)n 



h d(j)n 2/i I 



1 — ao 



(4.18) 



from where the above condition for AV follows readily. 
For the single solitary wave solutions of the sn. cn and 
dn type one wants to choose 



dAV ai(j)l + a24>l + as^r, 



(4.19) 



with fli dependent on the choice of elliptic function, in 
order to cancel the effect of the (/)^ terms in the corre- 
sponding equation of motion. 
Integrating, we get for AV 

AV =^ + ,^J_^f^ 



2a2 2(1 -a02) 

'0,1 02 

. a 



02 



a3\ 
a ) 



If we assume a soliton solution of the form 

0,1 = Aa sn|/3[(n + c)e — ut], ?ti| , 

where c is an arbitrary constant, then 

= l3^v^{Aa)^[l - (1 + m) sn^ + m sn^ 
which leads to 



Oi = — ma 



(iv 
Aa 



a-i 
as 



a/3^w^(H-m), 
- a{PvAa)'^. 



(4.20) 
(4.21) 
(4.22) 

(4.23) 

(4.24) 
(4.25) 



Explicitly, for the sn type of solitary wave with a kink 
limit when m = 1 we need to choose the extra potential 
term to satisfy 

a{l3vAa)^(j)n \^ {I + m)(f>l m4>l 



dAV 



1 



(Aa)2 (Aa)4 

(4.26) 

If instead, we assume a pulse soliton solution (for 
A < 0) of the form 



Aa dn 



|/?[(n + c)e-ut],m| , (4.27) 



then we have 



02 = l3'^v^{Aa)'^[{m - 1) + (2 - m)An^ - dn*], (4.28) 

which leads to the (different set) of coefficients {oj} for 
the extra term needed to be added to the Hamiltonian 



Qi 



a 



(3v' 
Aa, 

02 = — aP^v'^{2~m) 

03 = (1 — m)a{(3vAa)'^ 



(4.29) 

(4.30) 
(4.31) 



Explicitly, for the case of dn solitary waves we choose our 
additional potential term to satisfy: 



dAV a{PvAa)^(j)n 



d<t>n (1 - )2 



(4.32) 



(m- 1) + 



(2-m)02 04 
(Aa)2 (ylo)4 
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Finally, we assume a pulse solution (for A < 0) of the 
form 

(j)„ = Aa cn|/3[(n + c)e - + c], m| , (4.33) 

then we have 

= P^v^iAa)^ [{1-171) + {2m- l)cn^ -mcn^], (4.34) 

which leads to the (different set) of coefficients {a;} for 
the extra term needed to be added to the Hamiltonian 



(3v 



02 = - a/?\2(2m- 1) , 

03 = — (1 — m)a{Pv Aa)^ 



(4.35) 

(4.36) 
(4.37) 



Explicitly, for the case of cn solitary waves we choose our 
additional potential term to satisfy: 



d/W _a{f3vAaf 
d(j)„ 



(1 - «0D 
X (l-m) + 



(4.38) 

(2m - 1)^2 ^^4 



{Aay 



{Aa) 



Because of our "fine tuning" of AV, in all cases, the 
solitary waves effectively obey the second order differen- 
tial difference equation: 



1 — ad 



1+1 + - 2(?!)„) 



(4.39) 



As long as a is proportional to e^, we expect that this 
equation has the correct continuum limit. Next, we 
demonstrate this for the three cases explicitly. 



A. Positive A and kink solutions 



If we assume a solution of the form 



4>n = Aa sn|/3[(7i + c)e - wt], mj , (4.40) 

where c is an arbitrary constant we obtain matching 
terms proportional to [sn„+i + sn„_i] that 



sn2(/3e) = a 



(Aa)2 



(4.41) 



From the term linear in sn we then obtain 



\a^ = ^ [1 - dn(/3e)cn(/3e)] - (1 + m)f3^v^. (4.42) 

When w = 0, we get our previous static solution (3.16). 
This equation becomes in the continuum limit (e —>■ 0) 



Xa^ 



(1 + m)(l - u2) ' 



(4.43) 



Finally, matching the cubic terms in sn yields an equation 
for A namely 



AA^a^ = 2m 



sn2(/3e) 



- 2to/3\2 



From this we deduce that 

m[sn2(/3e)-/32^;2g2] 



A' 



[1 - dn(/3e)cn(/3e)] - (1 + 7n)e^P^v^ /2 ' 
In the continuum limit we get 

2m 



A' = 



1 + m 



(4.44) 



(4.45) 



(4.46) 



which agrees with our continuum result (3.17). If we 
divide (14.41(1 by (|4.44() we obtain: 



sn2(/3e) 



A 2[sn2(/3e)/e2_/32„2]- 

For small lattice spacing we get the result 

Ae2 



a 



2(1 -w2)' 



(4.47) 



(4.48) 



This shows that in order to have a solitary wave of ve- 
locity V, the discretization depends on 

When u = this reduces to the static case: 

aa^A^ = m sn^{Pe) , 
2m sn2(/3e) 



2 a2 



Xa^A 

Xa^ = ^ [1 - cn(e/3)dn(e/?)] , 
which implies that 

aa^ = 1 — cn(/?e)dn(/3e) , 
a = Xe^/2 , 

m sn2(/3e) 



A' 



1 - cn(/?e)dn(/3e) 



Localized mode: Let us now consider the m 1 limit 
in which case we get a localized kink solution 



where 



(pn = Aa tanh[/3(n -I- c)e — vt] 



A^l, aa^ = tanh (/3e) 
2 

Xa^ = — tanh^(/3e) - 20^ 



(4.49) 
(4.50) 



B. Negative A and pulse solutions 

If we assume a solution of the form 

(pn = Aa dn|/3([n + c]e - vt), m| , (4.51) 
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where c is an arbitrary constant we obtain matching 
terms proportional to (dn„+i + dn„„i) that 



sn2(/3e) 



-a{Aaf 



cn2(/3e) 

From the term hnear in dn we then obtain 



dn(/3e) 
cn2(/3e) 



- 1 



'2„,2 



- (2 - m)P^v 



(4.52) 



(4.53) 



When w = 0, we get our previous static solution (3.21). 
This equation becomes in the continuum limit (e 0) 



(2 -m)(l - w2)' 



(4.54) 



Finally, matching the cubic terms in dn yields an equa- 
tion for A namely 



-XA'a' 



2 sn2(/3e) 



cn2(/3e) 
From this we deduce that 



(4.55) 



A' = T 

dn(/3e) - cn2(/3e) - (2 - m)e2/32„2cn2(/3e)/2 

(4.56) 

In the continuum limit we get 



A' = 



(4.57) 



which agrees with our previous contimmm result (3.22). 
Dividing H4.52|l by (|4.55() we obtain 



a 
X 



e2sn2(/3e) 



sn2(/3e) - e2/32„2cn2(/3g) 



For small e we have 



Ae2 



2(1 -i;2)' 



(4.58) 



(4.59) 



When V ^ 0, a ~ Xe'^ /2. This exactly cancels the A(/)^ 
term in the equation of motion and we get the simple 
discretization for the time independent case. In the time 
dependent case wc again have the unusual result that the 
discretization needed is a function of the velocity of the 
solitary wave. 

If instead, we assume a solution of the form 



|/3[(ri + c)e-ut],m| , (4.60) 



Aa cn 



where c is an arbitrary constant wc obtain matching 
terms proportional to (cn„+i + cn„_i) that 



sn2(/?£) 



(4.61) 



From the term linear in cn we then obtain 
cn(/3e) 



^Xa' 



1 



{2m-l)P^v\ (4.62) 



When t> = 0, we get our previous static solution (3.25). 
This equation becomes in the continuum limit (e —^ 0) 



Aa2 



(2m - 1)(1 - -y2)' 



(4.63) 



Finally, matching the cubic terms in cn yields an equation 
for A namely 



^XA'a 



2„2 



2m sn2(/3e) 



2„,2 



£2 dn2(/3e) 
From this we deduce that 



A' 



2mP^v 



(4.64) 



m 


sn^iPe) - P'^v^e^dn^iPe) 




cn{Pe)-dn{Pe) 


- (2m - l)e2/32„2 


dn'(/3e)/2 



In the continuum limit we get 

2m 



A' = 



2m - 1 



(4.65) 
(4.66) 



which agrees with our previous continuum result (3.26). 
Dividing (|¥^ET|) by we obtain 



a 
X 



e2sn2(/3e) 



sn2(/3e)-e2/32w2dn2(/3e) 



(4.67) 



For small e wc have 



Ae2 



2(1 -i;2) 



(4.68) 



When V ^ 0, a = Xe'^/2. This exactly cancels the A0^ 
term in the equation of motion and we get the simple 
discretization for the time independent case. In the time 
dependent case we again have the unusual result that the 
discretization needed is a function of the velocity of the 
solitary wave. 

Localized mode: In the limit of m — > 1, both cn and dn 
solutions reduce to the localized pulse solution 

(h„ = Aa sech |/3[(7i + c)e- wt]| , (4.69) 



where 



-aa^A^ 
-Xa^ 



s\n\?{Pe) , 

4[cosh(/?e) - 1] -/?2v2 



(4.70) 



-A 0^24 sinh"(/3e)- 2/32^2 
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FIG. 5: (Color online) Scattering of kink-anti-kink waves. FIG. 6: (Color online) Scattering of dn pulse-anti-pulse 
For completeness, at t = 0, we also depict the time-derivative waves. For completeness, at t = 0, we also depict the time- 
of the initial wave function (in green) . derivative of the initial wave function (in green) . 



V. SCATTERING OF SOLITARY WAVES 

When wc have two soHtary waves colUding with op- 
posite velocity, then the equation of motion depends on 
whether we are considering sn, cn or dn type sohtary 



waves. In general we have 



1 



1 — M 



-4>n4>n - (1 - a0„) 



2(/)„] +A(/)„(a2 



<t>l) 



(5.1) 
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where the partial derivatives of Al/ are given by 
Eqs. iHSnil, and (g^HIl, for the three types of so- 

lutions, respectively. The AV term given by Eq. (|4.20|) 
adds an extra term to the energy conservation equation. 
The conserved energy is given by 



E : 



E 



fn+l 



2(1 

+ X{a^a - 1) 



2e2 

ln(l-a02) 



2a' 



^ a;- 

Za 
AV 



(5.2) 



Typical scenarios are depicted in Figs. [3 and IHl for 

the scattering case of two kink-anti-kink waves, and two 
pulse-anti-pulse waves, respectively. The solution of the 
kink-anti-kink "scattering" waves appears to correspond 
to a spatially localized, persistent time-periodic oscilla- 
tory bound-state (or a breather |23| ) with some ra- 
diation (phonons) at late times. Breathers are intrinsi- 
cally dynamic nonlinear excitations and can be viewed 
as a bound state of phonons. The scattering of pulse- 
anti-pulse is different: there is a flip after collision and 
relatively less radiation. 



VI. ENERGY OF SOLITARY WAVES AND THE 
PEIERLS-NABARRO BARRIER 

In a discrete lattice there is an energy cost associated 
with moving a localized mode by a half lattice constant, 
known as Peicrls-Nabarro (PN) barrier [isl [T9l | . In the 
m ^ 1 limit the elliptic functions become localized and 
become cither pulses, or kink-like solitary waves. We 
shall now show that, rather remarkably, the PN barrier 
vanishes in this model for both types of localized modes. 
In fact we prove an even stronger result, that the PN 
barrier is zero for all three periodic solutions (in terms of 
sn, cn and dn). If the folklore of nonzero PN barrier being 
indicative of non-integrability of the discrete nonlinear 
system is correct, then a la AL lattice, this could also 
be an integrable model. It would be worthwhile if one 
can prove this explicitly by finding the Lax pair for this 
model. 

We start from the conserved energy expression given 
by Eq. with AV of Eq. (|On|) . First of all, we shaU 
show that in the time-dependent as well as in the static 
cases, the conserved energy is given by 



B\n{l - a(t>l) + A 



(6.1) 



where the constants A and B vary for each case. We first 
note that in view of Eqs. (|4.1t)|) and H4.19|l . we have 



2(1 - a02 ) 



2a(l-a02) 



(6.2) 



Combining this term with AV of Eq. (|4.20|) in the energy 
expression (|5.2|l . we then find that the conserved energy 



is given by 



E 



A ai 



2a 



"2^ 



a 

A 
2^ 



IJn+lfPn 



(6.3) 



-3+^)ln(l-o02) 



Quite remarkably we find that in the case of all the 
(i.e. sn, cn as well as dn) solutions 



1 



A 

2^ 



fli 



= 0. 



(6.4) 



where use has been made of relevant equations in Sec. lIVI 
As a result the (j)^ term vanishes. Thus the conserved 
energy takes a rather simple form as given by Eq. 16. 1|) 
in the case of all of our solutions, and the constants A 
and B are given by 



A 

B 



2a 
Ao^ 
2a 



A 

2a2 



a2 
2a2 



(6.5) 
(6.6) 



Here 01,02,03 have different values for the three cases 
and are as given in Sec. IIVI In the special case of the 
static solutions, there is a further simplification in that 
the constant term also vanishes, since 01,2,3 are all zero 
in that case. 

Let us now calculate the PN barrier in the case of the 
three periodic solutions obtained in Sec. IIVI and show 
that it vanishes in all the three cases. Before we give the 
details, let us explain the key argument. If we look at 
the conserved energy expression as given by Eq. (|6.1|l . we 
find that there are two c-dependent sums involved here. 
We also observe that in these expressions, time t and the 
constant c always come together in the combination 



ki = (3{ce — vt) 



(6.7) 



Further, using the recently discovered identities for Ja- 
cobi elliptic functions we explicitly show that for all the 
solutions the first sum is c- independent. Since the total 
energy E as given by Eq. H6.1|) is conserved, its value 
must be independent of time t and since time t and the 
constant c always come together in the combination ki 
as given by Eq. H6.7(l , it then follows that the second sum 
must also be c-independent and thus there is no PN bar- 
rier for any of our periodic, and hence also the localized 
pulse or kink, solutions. 

In particular, the following three cyclic identities |23 
will allow us to explicitly perform the first sum in 
Eq. inn: 

TO sn(a:) sn(a:: + o)— ns(a) (6-8) 
X [Z{x + a)~ Z{x) - Z{a)] , 
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m cn{x) cn[x + a) =m cn(a) 

+ ds(a) [Z{x + a) 

and 



Zix) 



(6.9) 

Z{a)] , 



dn{x) An{x + a) =dn(a) (6.10) 
+ cs(a) [Z(a; + a) - Z(a:) - Z{a)\ . 



Here ^(x) = Z{x,m) is the Jacobi zeta function |21ll2/ 
In addition, we use the fact that 



N 

|z[/3e(n+l) + fci,m]-Z(n/3e + fci,m)| 



. 



.11) 



Let us now consider the sum as given by first term 
of Eq. (|6.1|) for the three cases. We first consider the 
kink-like solution which can also be written as 



0„ — Aa sn(n/3e + ki , m) . 
Using the identities given above, we obtain 



(6.12) 



N 

E 

n=l 



i4>n+l4>n) 



(Aa)2ns(/3e) 



N Z{Pe) . (6.13) 



For the dn pulse-like case we have instead 
0„ = Aa dn(n/3e + fci, m) , 
and the first sum in Eq. H6.1(l becomes 

N 



E 

n=l 



(6.14) 



(6.15) 



N{Aaf 



[dn(/3e) - cs(/3e)Z(/3e)] 



Finally, for the cn pulse-like case we have 
(j)n = Aa cn(n/3e + fci, 771) . 
In this case, the first sum in Eq. (|6.1|l becomes 



(6.16) 



E 



N{Aa) 



— Z{(3e)ds{(3e) - cn(/3e) 
m 



(6.17) 



It is worth noting that all these sums are independent 
of the constant fci and hence c. As argued above, since 
the total energy E as given by Eq. (|6.1() is conserved (and 
hence time-independent) and since t and c always appear 
together in the combination fci in Eq. (|6.1|1 therefore it 
follows that the second sum in Eq. H6.1|l must also be 
independent of fci and hence c. We thus have shown that 
the PN barrier is zero for the three periodic solutions and 
hence also for the localized solutions (which arc obtained 
from them in the limit tti = 1). 



In fact, it is easy to show that for all the three cases, 
the second sum (apart from a trivial fci-independent con- 
stant) in Eq. (6.1) is given by 



N 



El- 

71=1 



sn^(/3e) 
cn2(/3e) 



dn^(n/3e + fci) 



(6.18) 



In particular, for the solution H6.12|l the second sum is 
given by 



N 



bJ2\ ln[cn2(/3e)] +ln 



1 



sn2(/5e) 



dn^(n/3e + fci) 



cn2(/3e) 

(6.19) 

On the other hand, for the solution (|6.14|l . the second 
sum is given by 



B 



N 

Ei4^ 



sn2(/3e) 
cn2(/3e) 



dn^(n/3e -I- fci 



(6.20) 



Finally, for the solution H6.16|l . the second sum is given 
by 



5^ I In 

n=l 



cn2(/3e) 
dn2(/3e) 



In 



1 



cn2(/3e) 



dn^(ri/3e + fci) 



(6.21) 

It is worth remarking at this point that by following 
the above arguments, it is easily shown that even in the 
AL model the PN barrier is zero for the periodic dn 
and cn solutions. In particular, since in that case the 
energy is essentially given by the first sum in Eq. (|6.1|l . 
hence using Eqs. (|6.15|) and 16.17|) it follows that indeed 
for both dn and cn solutions jJl there is no PN barrier in 
the AL model. 

In general, we do not know how to write the sum in 
Eq. (|6.18|) in a closed form. However, for m ~ 1 and 
N —I- cx), the sum of logarithms in Eq. (|6.3|) can be carried 
out explicitly, and the energy of the solitary wave can be 
given in a closed form. 

We proceed as follows: One can show that, for N 
00, the following identity can be derived from the AL 
equation (see, for instance, Ref. (2^): 



r -1 2 

ln|l + sinh^(p) secli^[p(77-.T)]| = -. (6.22) 



Then, for kink-like solutions and m = 1, we have 

ln(l-a(/i^) ln[l - ianh\l3e) tanh^(n/3e + fci)] 

= lnsech^(/3e) + ln[l + sinh^(/3e) sech^r7/3e -t- fci)] . 

Thus, in the limit when N —^ 00 and 771 = 1, the sum of 
logarithms in Eq. 16.3|l becomes 



(3e 



J2 hisech2(/3e) 



(6.23) 
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Similarly, in the case of the two pulse solutions, dn and 
cn are identical for m = 1, and we can write 

\n{l-a(j)l) ln[l + smh^{(3e) sech^{n(3e + ki)] . 

In the limit when N ^ oo, the sum of logarithms in 
Eq. H6.3|) for the two pulse solutions, dn and cn, for 
TO ~ 1, is simply equal to 2/{(3e). 

For the case of the kink solution, for to = 1, B as given 
by Eq. simplifies to 

B = -^^^^, (6.24) 

whereas for the case of the pulse solutions, at to = 1, U 
simplifies to 

S = -2^^. (6.25) 

ae 

VII. CONCLUSIONS 

In this paper we have shown how to modify the naive 
discretization of field theory so that the discrete the- 
ory is a Hamiltonian dynamical system containing both 
static and moving solitary waves. We have found three 
different periodic elliptic solutions. For time-dependent 
solutions we found the unusual result that the discretiza- 
tion is dependent on the velocity. 

In the static case, we have studied the stability of both 
kink-like and pulse-like solutions, and have found differ- 
ent qualitative behavior of the lowest eigenvalue of the 
stability matrix in the two cases. For typical values of 
the model parameters, in the case of kink- like solutions, 
we found that stability requires the number of sites, 
to be larger than a minimum value, while for pulse-like 
solutions stability is achieved for arbitrary values of N. 
The magnitude of the lowest eigenvalues increases with N 
in the kink-like case, and decreases with N for pulse-like 
solutions. The lattice spacing, e, is not an independent 
parameter and always decreases with N. 



We also determined the energy of the solitary wave 
in the three cases. Using the Hamiltonian structure, we 
were able to argue that the PN barrier [TM Il9l | for all 
solitary waves is zero. This leads to the intriguing ques- 
tion whether this Hamiltonian system is integrable. As 
an additional result we explicitly showed that for the 
two elliptic solutions (dn and cn) Q of the integrable 
Ablowitz-Ladik model Q the PN barrier is zero-as one 
would expect. 

The single solitary wave solutions were found to be sta- 
ble and when we scattered two such single-kink solitary 
waves there were two different behaviors. For pulses, 
the pulse-anti-pulse solution leads to scattering with a 
fiip and a little radiation (phonons). The existence of 
phonons brings into question the integrability of this 
model, in spite of the zero PN barrier. For kink-like so- 
lutions we found a breather-like behavior [27l | during the 
collision. However, we have not found exact two-solitary 
wave or breather solutions, which would help clarify the 
integrable nature of this system. 

The results presented here are useful for structural 
phase transitions fl6l l2^ and possibly for certain field 
theoretic contexts [13 • Our results also hold promise 
for appropriate discretizations of other discrete nonlin- 
ear soliton-bearing equations P, Possible extension 
to discrete integrable models in 2 -I- 1 dimensions, e.g. 
Kadomtsev-Petviashvili (KP) hierarchies [s^l, would be 
especially desirable. Extension to time-discrete inte- 
grable models [sij is another interesting possibility. 
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